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A finite-size quasi two-dimensional Bose-Einstein condensate collapses if the attraction between 
atoms is sufficiently strong. Here we present a theory of collapse for condensates with the inter¬ 
atomic attraction and spin-orbit coupling. We consider two realizations of spin-orbit coupling: the 
axial Rashba coupling and balanced, effectively one-dimensional, Rashba-Dresselhaus one. In both 
cases spin-dependent “anomalous” velocity, proportional to the spin-orbit coupling strength, plays 
a crucial role. For the Rashba coupling, this velocity forms a centrifugal component in the density 
flux opposite to that arising due to the attraction between particles and prevents the collapse at 
a sufficiently strong coupling. For the balanced Rashba-Dresselhaus coupling, the spin-dependent 
velocity can spatially split the initial state in one dimension and form spin-projected wavepackets, 
reducing the total condensate density. Depending on the spin-orbit coupling strength, interatomic 
attraction, and the initial state, this splitting either prevents the collapse or modifies the collapse 
process. These results show that the collapse can be controlled by a spin-orbit coupling, thus, 
extending the domain of existence of condensates of attracting atoms. 

PACS numbers: 67.85.Fg, 67.85.Hj, 05.09+m 


I. INTRODUCTION 

Understanding Bose-Einstein condensates (BEG) of in¬ 
teracting particles is one of the most interesting prob¬ 
lems in condensed matter physics [l[ . For uniform three- 
dimensional systems repulsion between the bosons de¬ 
pletes the condensate, while attraction leads to the con¬ 
densate instability seen as the appearance of Bogoliubov 
modes with imaginary frequencies. For the finite-size 
condensates this instability can be seen in their collapse 
The collapse, where the size of the state goes to 
zero after a finite time, strongly depends on the spatial 
dimension D and is possible only in the D = 2 and D = 3 
condensates. The physics of the collapse is related to the 
fundamental problems of nonlinear optics and quantum 
mechanics Q, plasma instability and polaron forma¬ 
tion 0. 

The main features of the collapse of a free, not re¬ 
stricted by an external potential, condensate, are deter¬ 
mined by the interplay of its positive quantum kinetic 
and negative attraction energies dependent on the char¬ 
acteristic size of the condensate a. The kinetic energy 
is proportional to while the attraction contribution 
behaves as —a~^. For D = 3 the dependence of the to¬ 
tal energy on a is non monotonic and the collapse with 
a —>■ 0 occurs at any interaction strength since at small a 
the attraction dominates 0. For D = 2 the interaction 
and kinetic energies scale as and the collapse occurs 
only at a strong enough attraction. 

The BEG phy sics becomes much richer with synthetic 
gauge fi elds |ll| and synthetic spin-orbit coupling (SOC) 
[l2l Il3l|. For the latter, optically produced atomic pseu¬ 


dospin 1 /2 is coupled to atomic momentum and to a syn¬ 
thetic magnetic field. The SOG can be produced in var¬ 
ious forms, simulating the Rashba and the Dresselhaus 
symmetries Bin known in solid state physics. This 
coupling opens a venue to the appearance of new phases 
in a variety of ultracold bosonic [Tfil - I^ and fermionic 
[26l - l^ ensembles. It is well appreciated that the SOC 
plays crucial role in BEG physics in uniform three- 
dimensional gases with interparticle repulsion (Mm. 
For D = 2 the phases of the BEG of repelling bosons 
trapped in a harmonic potential were found in Ref. [s^ . 

One of the advantages of cold atomic gases is the fact 
that due to a very large particle wavelength compared 
with the atomic radius, the interatomic interaction can 
be accurately described by a single parameter, the scat¬ 
tering length Os, where positive (negative) Us corresponds 
to repulsion (attraction) between the atoms. The attrac¬ 
tion can be achieved by means of the Feshbach resonance 
[3^ in a certain range of the system parameters. Here 
we study joint effect of the interatomic attraction and 
spin-orbit on the spread and collapse of a quasi two- 
dimensional spin-orbit coupled BEG. 

This paper is organized as follows. In Sec. H we show 
by qualitative arguments, a variational approach, and di¬ 
rect numerical solution of the Gross-Pitaevskii equation, 
that the effect of the anomalous spin-dependent velocity 
due to the spin-orbit coupling [s^ can either completely 
prohibit the collapse or strongly modify the collapse pro¬ 
cess. We study the condensate dynamics and analyze the 
conditions at which the collapse does not occur. Possible 
relations to experiment and conclusions will be given in 
Sec. HI. 
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II. COLLAPSE IN THE PRESENCE OE 
SPIN-ORBIT COUPLING 


A. General formulation: Hamiltonian and the 
collapse process 


We consider a pancake-shaped condensate of pseu¬ 
dospin 1/2 particles described by a two-component wave 
function d' = {r, t), {r, t)] , where r = (x, y), nor¬ 
malized to the total number of particles ^ 1. In the 
presence of the spin-orbit coupling, the evolution of the 
wavefunction is described by a system of coupled nonlin¬ 
ear partial differential equations in the Gross-Pitaevskii- 
Schrodinger form 






- a) - g2\'Sf 


vk. (1) 


Here M is the particle mass, Hso is the SOC Hamiltonian, 
B is the effective magnetic field, and a = (a^x,^y,^z) is 
the spin operator. The coupling constant in Eq. m is 
given by g 2 = —^irfPPas/Maz^ which we assume for sim¬ 
plicity to be spin-independent, where is the condensate 
extension along the z— axis, and Og is negative 
Below we consider two strongly different forms of Hso'- 
the Rashba coupling with the spectrum axially symmet¬ 
ric in the momentum space and the balanced, essentially, 
one-dimensional Rashba-Dresselhaus coupling. 

Without loss of generality, we consider an initial state 
prepared in a parabolic potential at zero temperature as: 


'I'(r, t = 0) = A(0) exp 


2a2(0) 


^A(O), 


( 2 ) 


where ■0(0) is the initial spinor, H(0) = N/tt/ a{{))^ and 
a(0) is the initial width. At t = 0, the confining poten¬ 
tial is switched off [I| and the spin-orbit coupling and 
the attraction between the atoms are switched on. The 
subsequent dynamics is, thus, a response of the system 
to the instantaneous change in the potential, interaction, 
and spin-orbit coupling. 

In what follows we use the units H = M = 1 and the 
dimensionless interaction g 2 = —^iraa/az- The unit of 
length f can be chosen arbitrarily, and the corresponding 
unit of time is 

We address first the collapse without spin-dependent 
effects. Here the energy of the system is 


E = -- 


-I- 52 l^'l 


dxdy, 


(3) 


and the evolution can be described by a variational ap¬ 
proach based on Gaussian ansatz Q 


'I'(r, t) = A{t) exp 



(1 -k iK{t)) 


i’i.O), 


(4) 


where the variational parameters b^{t) and a^{t) are the 
chirp and the packet width, respectively. The equa¬ 
tion of motion for Ov becomes Ov = — A/a^, where 


o 
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FIG. 1: (Color online) Time dependence of the condensate 
width for 52 N = IGtt and the values of a marked near the 
lines. The green short-dashed line corresponds to the absence 
of collapse. 


A = {g 2 Ef — Av) /2. The collapse occurs if exceeds 
the variational threshold value Av = Stt [^. The solu¬ 
tion of this equation is 


I At2 

av(t)=«(0)Wl-^. (5) 

The timescale of the evolution is the collapse time Tc = 
a^(0)/-\/A, and the characteristic collapse velocity is Vc = 
a(0)/T, = VX/a{0). 

The key point in the understanding of the role of the 
spin-orbit coupling in the collapse process is the modified 
velocity 

v = kH-VkHso, (6) 

with k = —id/dr, including the anomalous velocity 
term Vki?so (here Vk = d/dk) directly related to the 
particle spin. The evolution of the probability density 
p = ig given by the continuity equation 

|^ + V-J(r,t)=0, (7) 

with the components of the flux density 


J(r,t) = - 


Vkffs, 


T-. (8) 


The spin components of the condensate are given by ex¬ 
pectation values 


(Mt)) 


1 

N 


'k^di'i'dxdy. 


(9) 


B. Rashba spin-orbit coupling 

The first form of spin-orbit interaction that we consider 
is the Rashba coupling 

Hso = H]^ = (X {kx^y ky(7x) , 


( 10 ) 
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with the coupling constant a and k = ky). The corre¬ 
sponding spin-dependent terms in the velocity operators 
in Eq. m become 


dHR ^ OHr 
dk^ dky 


OlGx ■ 


( 11 ) 


The spatial scale of the SOC effects is described by the 
characteristic distance the particle has to move to flip 
the spin, = l/a. The corresponding spin rotation 
angle at the particle displacement L is of the order of 
L/Lso- At the initial stage of the BEC evolution t <^Tc 
we obtain from Eq. ([T]) for = 0) in Eq. ([2|) with 

t^(o) = [i,oF, 


2a2(0) ■ 

( 12 ) 

As a result, the component begins to grow at 

distances r ^ a(0) with a rate proportional to a. At a 
sufficiently large a this growth can eventually lead to the 
collapse prevention. 

At t > 0, the spatially nonuniform spin evolution be¬ 
gins. Since the spin precession angle at the displacement 
of a(0) is of the order of a(0)/Tso, starting from the fully 
polarized ipiO) = [1,0]^ state, the atoms acquire the 
anomalous velocity of the order of a(0)/Lso x a ^ Q!^a(0) 
for the weak SOC a(0) ^ Lso, or of the order of a other¬ 
wise. The criterion of a large spin rotation in the collapse 
is a(0) > Lso, that is a > l/a(0), while the condition of a 
sufficiently large developed anomalous velocity is a > z;c, 
that is a > \/A/a(0). If the latter inequality is satisfied, 
the centrifugal component in the flow caused by the SOC 
can prevent the collapse, as we explain in detail be¬ 
low. The condition of a weak effect of magnetic field on 
the collapse can be formulated as smallness of spin pre¬ 
cession angle due to the Zeeman splitting compared to 
the precession angle due to the spin-orbit coupling, that 
is TcB "C min {a(0)/Lso, 1} ■ We will assume this condi¬ 
tion and neglect the effects of the Zeeman splitting. 

We begin the analysis of the joint effect of the SOC 
and the interatomic attraction with numerical results ob¬ 
tained by direct integration of Eq. © for a strong attrac¬ 
tion, g 2 N ^ 1, where the effect is clearly seen, taking 
the initial spin state ■0(0) = [1,0]^. Eigure [T] shows the 
time-dependent width of the packet defined as 



where T is obtained by a direct solution of Eq. ([T]) for 
several values of a. The solid line in Fig. [T] corresponds 
to the collapse at a = 0 where in the vicinity of Tc, the 
numerically calculated using Eqs. m and (USD, width 
a{t) is accurately described by variational Eq. ([5]) with 
alt) (Tc - 

When spin-orbit coupling is included, the following 
features may be seen, (i) At short time t <C Tc, the 




0'*'(r, t —>• 0) = f 




x + iy 


V^a3(0) Ts, 


■ exp 



r/a(0) 



FIG. 2: (Color online) Density profile p{r,t) for 5^2A = IGrr. 
Black solid line is for t = 0, red dashed line is for t — 0.2a^(0), 
blue dot-dashed line is for t = 0.4a^(0), and green dotted 
line is for t = a^(0). (a) Here a — 0.67uc, and the dotted 
green line shown in detail in the inset, clearly demonstrates 
the collapse, (b) Here a = O.Sdric, the condensate is robust 
against attraction and can spread without collapsing. 


attraction-induced velocity develops linearly with t, while 
the anomalous velocity increases as As a result, the 
a (f)—dependences for all values of a are the same at small 
t. (ii) The packet width a{t) increases with time, reaches 
a plateau, and then decreases to zero. Thus, with the 
increase in a, the collapse still can occur, albeit taking 
a longer actual time tc > Tc- (iii) Increasing further, a 
reaches a critical value acr ~ 0.7vc such that at a > Ocr 
the anomalous velocity is large enough to prevent the 
collapse. The dependence of tc on the SOC strength can 
be described as tc ~ (ocr — cr) 

To get an insight of the effects of SOC on the collapse, 
we depict the density profiles in Fig. [21 At a large a 
the density forms a double peak with the maxima posi¬ 
tions separating with time as a result of the centrifugal 
component in the flux. The resulting two-dimensional 
density distribution is given by a ring of radius R(t) and 
width w(t) with a(t) ^ y/R{t)w(t), responsible for the 
broad plateaus in a{t)/a{0) seen in Fig. [T]at subcritical 
spin-orbit coupling. At R{t) ^ a{t) the interatomic in¬ 
teraction energy tends to zero as —1/R{t)w{t), and the 
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conserved total energy is the sum of the kinetic and SOC 
terms. At a < (see Fig. UKa)) the attraction is 
still strong enough to reverse the splitting and to restore 
the collapse. At a > acr, (see Fig. HKb)) the anoma¬ 
lous velocity takes over, the splitting continues, and the 
collapse does not occur . This process is naturally ac¬ 
companied by evolution of the condensate flux and spin 
presented in the Supplemental material [s^ . 

C. Balanced Rashba and Dresselhaus couplings 

In this subsection we consider a one-dimensional cou¬ 
pling 



— ^RD — 


(14) 


which is equivalent to the balanced Rashba and Dressel¬ 
haus contributions and gauged out by an a;—dependent 
spin rotation U = exp [idzx/Lso\- 

For simplicity we consider the initial state correspond¬ 
ing to the spin oriented along the x—axis with ■0(0) = 
[1, /^/2. Due to the anomalous velocity dHno/dkx = 

OLOz (cf. Eq. (HU), the initial state splits into two spin- 
projected wavepackets moving in the absence of interac¬ 
tions with velocities ±a. As a result, the effective interac¬ 
tion decreases, and the collapse can be prohibited by this 
decrease. This happens, however, only at certain condi¬ 
tions, which we establish here. For qualitative analysis 
we use the ansatz 


'(/’'''’■‘•(rzp,^) = A(t)exp 



^icv{t)x , 


(15) 

where the upper (lower) sign corresponds to spin up 
(down) and position rip = (x =F dv{t),y), and in addi¬ 
tion to the variational chirp bv{t) and width we 

introduced the variational momentum Cv(t). From this 
ansatz we obtain, using an approach similar to that of 
Ref. [5, equations of motion for and dv- 


d-y 


92^ dy 
TT d^ 


exp 



(16) 


92 N 


1 + 1 - 


2dl 


exp 


2dl 


for given av(0) = a(0) and other initial conditions 


dv(0)=0, dy{0) = a, Sv(0) = 0, (17) 


where dv{0) is due to the anomalous velocity term lead¬ 
ing to the spin-dependent splitting. These equations 
show that the collapse disappears if the coupling is strong 
enough to sufficiently separate the spin components, that 
is at a certain time dy changes sign from negative to pos¬ 
itive. 


FIG. 3: (Color online) Time dependence of the condensate 
width for g 2 N = Stt in the presence of balanced Rashba- 
Dresselhaus SOC and different values of a as marked near 
the lines. 


Qualitative conditions of the collapse in the presence of 
spin-orbit coupling in Eq. (HI, which can be found from 
Eq. (ITSl) . are as follows. ltg 2 N > Air, the collapse always 
occurs since even if the spin states are well-separated, 
each of them still has the sufficient number of atoms. 
Depending on the interatomic interaction and SOC, one 
can either obtain the collapse at the origin, producing 
a spin non-polarized condensate, or two spatially sym¬ 
metric ones producing z—axis polarized condensates. If 
g 2 N < Att, the collapse occurrence depends on the SOC 
strength. 

At a sufhciently strong SOC, the spin splitting of the 
initial state and possible collapse happen on different 
time scales. The splitting occurs fast, on the timescale 
of a(0)/a, and the interatomic attraction starts to play a 
role after the splitting. The condition of time scale sep¬ 
aration, which allows one to treat the splitting and the 
collapse independently, is formulated as a(0) < Tea or, in 
other words, as a > •\/A/a(0). This looks similar to the 
above condition for the critical Rashba coupling. How¬ 
ever, these conditions are qualitatively different. For the 
Rashba coupling, the density decreases to zero and the 
collapse disappears completely at any SOC stronger than 
the critical one. For the balanced Rashba-Dresselhaus 
coupling the maximum density decreases at most by a 
factor of two, and, therefore, the collapse can occur even 
at a very strong SOC, when spin-up and spin-down states 
are already well-separated in space. 

Figure [3] shows the time dependence of the packet 
width in Eq. m obtained by solution of Eq. dU) 
with spin-orbit coupling Hamiltonian (11411 for g 2 N = Stt. 
The behavior at small t Tc here depends on a since 
the peaks in the spin-projected densities split by 2a t 
due to the anomalous velocity. The numerically ob¬ 
tained critical value of a here is approximately 0.83xc 
and a(t) ^ {tc — t) shows a linear rather than a square- 
root behavior near the collapse time. 
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III. RELATION TO EXPERIMENT AND 
CONCLUSIONS 

To make connections to possible BEC experiments, we 
return to the physical units and estimate the constant 
g 2 as 0.05 for —as ^ lOOas ^ 5 x 10“^ ^m and ^ 1 
fim. The condition g 2 N > 2tt can be satisfied already 
for ensembles with N ~ 100 particles. The velocity of 
the collapse is Vc ~ hy/g^/Ma{Q). At a(0) ~ 10 ^m 
and N ~ 10^ this estimate yields Vc ~ 0.03 cm/s and the 
corresponding time scale Tf. = a(0)/uc ~ 0.3 s. Such a 
small value of Vc demonstrates that even a relatively weak 
experimentally achievable coupling can prevent the 
BEC from collapsing. At these conditions, the character- 

istic distance between the particles (^a'^{0)az/N^ ~ 0.5 

fim is much larger than —dTrOg < 0.1 /rm, still preventing 
a strong depletion of the condensate. 

To conclude, we have demonstrated that the anoma¬ 
lous spin-dependent velocity determined by the spin- 
orbit coupling strength can prevent collapse of a nonuni¬ 
form quasi two-dimensional BEC [4l|,|42. Eor the Rashba 
coupling with the spectrum axially symmetric in the mo¬ 
mentum space, this velocity leads to a centrifugal com¬ 
ponent in the two-dimensional density flux. As a re¬ 
sult, spin-orbit coupling can prevent collapse of the two- 
dimensional BEC if this flux is sufficiently strong to over¬ 
come the effect of interatomic attraction. In this case, the 
attraction between the bosons cannot squeeze the initial 
wavepacket and force it to collapse. In the case of ef¬ 
fectively one-dimensional balanced Rashba-Dresselhaus 
couplings, the anomalous velocity splits the initial state 
into spin-polarized wave packets, decreases the conden¬ 
sate density and, thus, can prevent the collapse. Our 
approach can be generalized in a straightforward way for 
intermediate case, where Dresselhaus and Rashba cou¬ 
plings have different strength. These results show that 
one can gain a control over the BEC collapse process 
by using the experimentally available synthetic spin-orbit 
coupling fields and, thus, extend the experimental abil¬ 
ities to study various nontrivial dynamical regimes in 
Bose-Einstein condensates of attracting particles. 
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Supplemental Material for “Collapse of spin-orbit coupled Bose-Einstein condensate” 


V. GRAPHIC MATERIAL FOR STRONG INTERATOMIC INTERACTION. 


Here we present additional graphic material to illustrate the role of spin-orbit coupling in the physics of the collapse 
of two-dimensional Bose-Einstein condensates. 
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FIG. 4: Flux of the condensate, defined in Eq. (8) of the main text, at 5^2^ = Idvr, a — 0.84nc, and t = 0.2a^(0). Since here 
a > Ocr = 0.7vc, this is the no-collapse regime, corresponding to a plot in Fig. 2(b) of the main text. This flux distribution 
leads to the increase in the width a{t) as defined in Eq. (13) of the main text. 



t/aHO) 


FIG. 5: Time dependence of the total condensate spin component as defined by Eq. (9) in the main text. Here g 2 N = Ifiyr, 
black solid line is for a = 0.67vc < eta (collapse regime), and red dashed line is for a = 0.84nc (no-collapse regime). 
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VI. NEAR-THE-THRESHOLD COLLAPSE 


Now we address a near-the-threshold collapse, which takes a long time, and where the difference between the 
variational g 2 N = 2tt and the exact = g 2 N = 1.8627r threshold couplings becomes important. We take g 2 N = 27r, 
where in the absence of spin-related effects, the total energy in Eq. (3) of the main text of a Gaussian state is 
zero. In this case the critical spin-orbit coupling sufficient to destroy the collapse by the anomalous velocity is 
determined by condition ~ Vc/a{0) and can be estimated as {g 2 N — Aex)^'^^ • The numerically obtained critical 
value is Ocr ~ 0.38(27r — Aex)^^^- Figure [5] shows that for a slightly larger than the critical value, the condensate 
first narrows and then broadens. The peak structure of Fig. 2 in the main text is not formed here, and the collapse 
disappears due to the broadening rather than due to the splitting. Although the change in the total spin shown in 
Fig. [7] is moderate compared to that presented in Fig. [SJ as expected for relatively small values of a, the collapse 
does not occur here. 



x/a{0) 


FIG. 6: (Color online) Density profile of the condensate for different times, 5^277 = 27r; a = 0.39(27r — Aex)^^^ (slightly above 
the critical value), black solid line is for t = 0, red dashed line is for t — a^(0), and blue dot-dashed line is for t = 4.4a^(0). 
Here no collapse occurs. 



FIG. 7: Time dependence of the total condensate spin component as defined by Eq. (9) in the main text. Here g 2 N = 27r, 
black solid line is for a = 0.19(27r — Aex)^^"* (collapse regime) and red dashed line is for a = 0.39(27r — Aex)^^^ (no-collapse 
regime). 
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